Manifestation of the Nuclear Anapole Moment in Ml Transitions in Thallium 
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We calculate nuclear spin-dependent parity non-conserving .El-amplitudes for optical transition 
6pi/2,F — * 6p3/2,F' an d for hyperfine transition &P\/2,f —* 6pi/2,F' m 205 T1. Experimental limit 
on the former amplitude placed by Vetter et al. [PRL,74,2658 (1995)] corresponds to the anapole 
moment constant n a = —0.26 ± 0.27. 
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In 1980 Flambaum and Khriplovich JlJ pointed out 
that the nuclear spin-dependent (NSD) part of the par- 
ity non-conserving (PNC) interaction in heavy atoms is 
dominated by the contribution of the anapole moment 
(AM) of the nucleus [Q] . After that AM was observed in 
the PNC experiment with 133 Cs ||, where the measured 
value of the AM constant n a appeared to be even larger 
than theoretical prediction for the "best values" of the 
constants of the nuclear PNC interaction (see Q and ref- 
erences therein). On the other hand, in the most accurate 
measurement of the PNC amplitudes 6pi/2,F ~ * &P3/2.F' 
in 205 Tl the NSD amplitude was found to be consis- 
tent with zero and smaller than theoretical predictions 

Si- 

In Ref. H the ratio between NSD amplitude and the 
dominant nuclear spin-independent (NSI) PNC ampli- 
tude was calculated in the one-particle approximation. 
Here we recalculate this ratio using CI+MBPT method 
||, [|, which allows to account for both core- valence 
and valence-valence correlations. We found that correla- 
tion corrections are relatively large, but do not explain 
the discrepancy between the measurement and the 
theory [Q. A more accurate measurement of NSD am- 
plitude in the optical transition 6p 1 / 2 ,F ~ * §Pz/2,F' is 
hampered by the much larger NSI amplitude and by the 
smallness of the hyperfine structure of the upper state. 
Consequently, it may be easier to measure PNC am- 
plitude for the hyperfine transition &P\/2,f ~~ * 6pi/2,F'> 
where NSI amplitude turns to zero, while NSD amplitude 
is not suppressed JTcfl . Here we find correlation correc- 
tions to this amplitude to be 20%. 

In the PNC experiments on the 6pi/2 - * 6^3/2 transi- 
tion in Tl, the ratio: 



K = Im(£l PNC /Ml) 



(1) 



of the PNC amplitude to the magnetic amplitude was 
measured with 1% accuracy in Ref. D and with 3% ac- 
curacy in Ref. pljl . In those experiments the hyperfine 
structure of the lower level Qpi/2 was resolved. That al- 
lowed to determine H(F) for two hyperfine levels F = 
and F = 1 of the ground state. For the level F = the 



only transition to F' = 1 of the level 6^3/2 is allowed, 
while for the level F = 1 transitions to both upper hy- 
perfine levels F' = 1,2 are allowed. Accordingly, 71(1) is 
some average for two transitions: 

11(0) = 71(0,1), (2a) 
11(1) = x 2 1Z(1,1) + (1-x 2 )1Z(1,2), (2b) 

where coefficient x 2 depends on the intensity of the tran- 
sitions and on experimental conditions, such as line width 
and optical depth. 

Observation of the independence of the PNC ampli- 
tude is important as it can give information about NSD 
part of the PNC interaction: 



H 



PNC 



\si +H NSD = "^|(--y-75 + ~dI)/AF). !>!) 



where Gf = 2.2225 x 10~ 14 a.u. is the Fermi constant of 
the weak interaction, Qw is the nuclear weak charge, k is 
the dimensionless coupling constant, 75 and a = 707 are 
the Dirac matrices, / is the nuclear spin (I — i for both 
stable isotopes 205 T1 and 203 T1), and p(r) is the nuclear 
density distribution. 

There are three main contributions to the coupling 
constant n in NSD part of the PNC interaction (|): 



«2 + kq v 



(4) 



where AM contribution is given by the constant n a [[[), 
the constant K2 — ^(4 sin 2 6\y — 1) ~ —0.06 corresponds 
to the NSD weak neutral currents f2~ij| . The term kq w is 
induced by the interference of the NSI interaction with 
the hyperfine interaction. For heavy nuclei this constant 
is proportional to A 2 / 3 @ ||, and for Tl k Qw w 0.02. 
Substituting these numbers in (0), we get: 



-§(*, 



0.06) . 



(5) 
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Theoretical predictions for AM constant depend on nu- 
clear model and vary within the range 0.1 < K a < 0.4 
(see H and references therein). On the other hand, for 
a given nuclear model, one can use measured values of 
n a to get information on the coupling constants of the 
nuclear P-odd interaction |lj, |lj, 1(| . 
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In this article we calculate the NSD amplitudes 
6pi/2,F ~~ * 6p 3 / 2 ,_F' and use Eq. (||) and experimental re- 
sults fromj5) to place a limit on the AM constant n a . 
Following |^7^| we use parametrization: 



K{F,F') = C(Z)[Q W - 6k£(F,F% 



(G) 



that links NSD amplitude to NSI one via the function 
£(F, F'). According to (||), one can define a function 
£(F) as follows: 

m = £(0,1), (7a) 
£(1) = z 2 £(l,l) + (l-z 2 )£(l,2). (7b) 

An important property of the one-particle approxima- 
tion is the equality £(1,1) = £(1,2) ||, which means 
that £(1) does not depend on the coefficient x 2 in ( pb| ) 
and (7b). Numerical values, obtained in (6), are: 



£ op (0) = 0.87; £ op (l) = -0.29. 



(8) 



In general, when electron correlations are taken into ac- 
count, £(1, 1) 7^ £(1, 2). Then, one has to use Eq. (^J) for 
TZ(F 7 F') and calculate function £(F', F). After that, ex- 
perimental function £(F) is given by (0). Consequently 
the separation of NSI and NSD amplitudes depends on 
the factor x 2 . 

NSI amplitude was studied many times, the most ad- 
vanced and accurate calculations being jl7], [l8| (for ear- 
lier references see |n|). It was shown there that many- 
body corrections to PNC amplitudes in Tl can be impor- 
tant. That stimulated us to recalculate function £(F, F'). 
We follow here the same procedure, which was used 
in iQ. It is based on the combination of the many- 
body perturbation theory for core-valence correlations 
and the configuration interaction for three valence elec- 
trons (CI+MBPT method) @,|, §. 

Most of the technical details of this calculation, such 
as basis sets, configuration sets, etc., are the same as in 
Ref. |Q , where a number of test calculations were made 
for the spectrum, hyperfme constants, El-amplitudes, 
and polarizabilities. All these parameters were shown to 
be in good agreement with the experiment. That allowed 
to estimate the accuracy of the calculation of NSI ampli- 
tude to be better than 3%. Here we use the same wave 
functions for the states 6pj, but neglect several smaller 
corrections, such as structural radiation, to the effective 
operators for valence electrons. The normalization cor- 
rection is the same for NDI and NSD amplitudes and 
does not affect the function £(F, F 1 ). 

In order to find PNC amplitude we solve inhomoge- 
neous equations: 



(E 6p3/2 - H eS )*i D l = Df * 6pi 



/2, 



(E 6pi/2 - ff ett )* 



6p 3 /2, 



(9) 
(10) 



where H is the effective Hamiltonian for valence elec- 
trons, which accounts for core- valence correlations within 
the second order many-body perturbation theory [Q, ^J, 



_D° is z-component of the effective El-amplitude in the 
length-gauge ^0|, and m is magnetic quantum number. 
Solutions of these equations can be decomposed in terms 
with definite angular quantum number J: 



(D) 



{D) ■ i-ab 



(11) 



NSI amplitude can be found by calculating the follow- 
ing matrix elements: 



1 



EInsi = (-l) 5 ~ m 2 n 2 ( 12 ) 



where we skip index m in matrix elements and take ad- 
vantage of the fact that -H^si * s diagonal in quantum 
number J. NSD part of the PNC interaction (||) can 
change this quantum number, and corresponding ampli- 
tudes have more complicated form: 



El 



NSD 



5/2 

£ 

,7=1/2 



C(J,F,F>)U* 6p3 JHf SD \¥ a D }) 



+ (^SWsdI^ 



(13) 



where constants C(J, F, F') are some combinations of the 
6j-coefficients (see |2l| for details). 

All wave functions in Eqs. (12), ( |l3| ) are many-electron 
ones. In the one-particle approximation, these expres- 
sions are simplified, and both NSI and NSD parts of the 
PNC amplitude have the form: 



El 



pnc — 



y-' (6p 3 /2 \\D\\ns 1/2 ) (ns 1/2 |g PNC |6jji/ 2 ) ^ , ^ 



£ ns 1/2 



The sum here runs over occupied (n = 1, . . . , 6) and va- 
cant (n > 6) states. Contribution of the occupied states 
with n < 5 is very small, while n = 6 contributes almost 
as much as the whole sum over vacant states. The term 
n = 6 corresponds to amplitudes with the index b in (12) 
and (fl3"l). It is seen, that all intermediate states in (14) 
have J = 1/2. This leads to the equality £(1, 1) = £(1,2), 
which is not correct for a more general case of Eq. (|l|) . 
The many-body corrections are strongest for the weak 
amplitude F = 1 — > F' = 1, which affects the value of 

£(!,!)■ 

Our results for the function £(F, F') are given in Ta- 
ble I. We find them from the calculated amplitudes jl^ ) 
and (|l5| ) in two approximations. At first, we use con- 
figuration interaction method for three valence electrons 
with conventional operators. Then, we use second or- 
der many-body perturbation theory to construct effective 
Hamiltonian H cS and random phase approximation for 
the effective operators D° ff and Hp^ c . 

It follows from the comparison of Table 1 with the 
one-particle approximation (H), that correlation effects 
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TABLE I: Table 1. Calculated values of £(F,F') in different 
approximations: configuration interaction (CI) for three va- 
lence electrons, and CI+MBPT method; a and b correspond 
to two contributions in Eqs.@ and 







CI 




CI+MBPT 


F, F' 


a 


b 


total 


a b total 


0, 1 


1.09 


1.29 


1.20 


1.08 1.12 1.10 


1, 1 


-0.498 


-0.513 


-0.506 


-0.500 -0.431 -0.462 


1, 2 


-0.337 


-0.413 


-0.378 


-0.331 -0.361 -0.348 



enhance NSD amplitudes. For the weakest amplitude 
F = 1 — ► F' = 1 the correlation correction exceeds 50%. 
For two other amplitudes correlations are less important, 
but still account for 20% - 25% enhancement. Valence 
correlations are larger for the amplitudes b. The domi- 
nant contribution to these amplitudes corresponds to the 
intermediate states from the configuration 6s6p 2 , where 
correlations between two p-electrons are very strong. In 
contrast to that, the main contributions to amplitudes 
a correspond to configurations 6s 2 np, where correlations 
are much weaker. 

We showed above, that correlation corrections to NSD 
amplitudes are rather large. Moreover, our values of 
£(1, 1) and £(1, 2) noticeably differ from each other. That 
leads to the dependence of the experimentally observed 
amplitude (^) on x 2 . The value of this parameter de- 
pends on the experimental conditions. In the linear 
regime, x 2 and 1 — x 2 are proportional to the intensi- 
ties of the corresponding lines. That gives x 2 = i [ Bj] . 
Actual experiment ]5| was done in the nonlinear regime, 
when in the center of the line the light was completely ab- 
sorbed, and PNC signal was detected only on the wings. 



In these conditions, one can expect that \ < x 2 < 
Below, we will perform analysis for each of the limiting 
cases. 

If we substitute values from Table 1 to Eq. (f7t]), we 
get: 



£(0) = 1.10; £(1) 



-0.367, x 2 = 1/6, r . 
-0.405, x 2 = l/2. (Lb) 



NSI amplitude can be found as weighted average: 
£(Q)tt(l) - ((l)K(Q) 



£(0)-£(l) 



(16) 



Experimental difference between 7Z(1) and 71(0) is only 
about 1%. Because of that, both values of £(1) from 
(H) lead to the same value of IZnsi = —14.68 x 10 in 
agreement with the result from || . 

The difference ATZ = 1Z(1) — 71(0) can be written as: 

aK . fc «a_«i! RNSI (17) 

ww 

= -4(k„ + 0.06) ^ (0) ~ ^ (1) ^nsi, (18) 
Ww 

where we use relation (||) between k and n a . Table 1 and 
Eq. (|lj) give £(0) - i(l) = 1.49 ± 0.02, and substituting 



the standard model value Qw = —116.7 |23[ , we get: 

ATI = (0.051 ± 0.001) (« a + 0.06)ft NS i, (19) 

where the error bar corresponds to two values of x 2 
in ([[5]) and does not account for the theoretical error, 
caused by the neglect of the higher orders of the many- 
body perturbation theory. The latter was estimated in 
ILsf for NSI amplitude to be close to 3%. Here we neglect 
the structural radiation corrections and few other correc- 
tions, which can contribute on the percent level, so we 
estimate the actual accuracy of Eq. ( [l9| ) to be about 5%. 
On this level, the uncertainty in experimental conditions 
described by the parameter x 2 is negligible. 
Using experimental values from [p|: 



ATI 



(-14.68 ±0.06 ±0.16) x 10~ 8 , 
(0.15 ±0.13 ±0.15) x 10~ 8 , 



we get the following result for the AM constant: 



-0.26 ±0.27. 



(20) 
(21) 



(22) 



In an independent measurement |Tlj of PNC effects in 
Tl a very close central value for the parameter ATZ was 
obtained, though with a three times larger uncertainty. 
If we use (§) instead of (|l|), we get n a = -0.32 ± 0.35. 
It means, that correlations account for 30% corrections 
and lead to a smaller absolute value of the AM constant. 
Note, that in Ref. Q the approximate values £(0) = 1 
and £(1) = were used instead of the more accurate 
one-particle values (i), and the relation k = — |k q was 
used instead of Eq. ( 5 ) . 

It was first recognized by Novikov and Khriplovich flio| , 
that NSD operator also leads to the .El-amplitude be- 
tween hyperfine sublevels of the same electronic state. 
The most interesting in this respect is the hyperfine tran- 
sition in the ground state. Such amplitudes were calcu- 
lated in the one-particle approximation for Cs and Tl 
JTo| and for K |^2| . The only many-body calculation was 
done recently for Fr pl[ . It is straitforward to rewrite 
Eq. (|l|) for this case, and all calculations are similar to 
those for optical transition. The result in a.u. is: 

(6p 1/2 ,l||M NSD ||6p 1/2 ,0) =2.11x 10- n i K , (23) 

where we use the same level of approximation, as above 
and add normalization correction |18| . In the one-particle 
approximation the Ml-amplitude for this transition is 
equal to — Correlations change this value only at 
sub-percent level, and we can safely use it to calculate 
71: 



^hf(6pi/ 2 ) 



-1.00x10~ 8 k (24) 
0.67 x 10~ 8 (k q ±0.06). 



Comparison of this value with the one obtained in |lC| ] 
shows that correlations increase the answer by approx- 
imately 20%. Result (p4|) can be compared also to the 
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F = 4 — > F' = 5 transition in the ground state 7s of 
211 Fr, where K = 3.0 x 1CT 9 k |§l). Though the Ml- 
amplitude for the hyperfine transition in Tl is signifi- 
cantly weaker than in Fr, it should be much easier to 
do the experiment with stable Tl, than with radioactive 
Fr. Note, that for lighter Cs, 1Z is an order of magnitude 
smaller. 

We see that electron correlations account for substan- 
tial corrections to AM amplitudes, but do not explain 
the difference between the experiment ^ and the nuclear 
theory prediction that n a = 0.25 ± 0.15 Q. The exper- 
imental accuracy for NSD amplitude for 6pi/2 —> 6P3/2 
transition is not high, because this amplitude is much 



weaker, than NS1 one. Therefore, it may be very inter- 
esting to measure the hyperfine amplitude (|23|), where 
PNC effects are completely determined by the NSD part 
of the weak interaction. Note also, that the frequencies of 
the hyperfine transitions for two natural isotopes 203 T1 
and 205 T1 differ by 1% and should be easily resolved. 
That gives the possibility to measure AM constants for 
each of the isotopes. 



The author is grateful to V. V. Flambaum, I. B. Khrip- 
lovich, and S. G. Porsev for valuable discussions. This 
work was supported by RFBR, grant No 02 - 02 - 16387. 



[1] 

[2] 
[3] 



[4] 

[5] 

[6] 
[7] 

[8] 

[9] 

[10] 

[11] 
[12] 



V. V. Flambaum and I. B. Khriplovich, Sov. Phys.-JETP 
52, 835 (1980). 

Y. B. Zel'dovich, Sov. Phys.-JETP 6, 1184 (1957). 
C. S. Wood, S. C. Bennett, D. Cho, B. P. Masterson, 
J. L. Roberts, C. E. Tanner, and C. E. Wieman, Science 
275, 1759 (1997). 

V. F. Dmitriev and I. B. Khriplovich, On the t heory 



of the nucl ear anapole moments (2002), E-print; nucl- 



th 020104] 



P. A. Vetter, D. M. Meekhof, P. M. Magumder, S. K. 
Lamoreaux, and E. N. Fortson, Phys. Rev. Lett. 74, 2658 
(1995). 

I. B. Khriplovich, Phys. Lett. A 197, 316 (1995). 

V. A. Dzuba, V. V. Flambaum, and M. G. Kozlov, JETP 

Lett. 63, 882 (1996). 

V. A. Dzuba, V. V. Flambaum, and M. G. Kozlov, Phys. 
Rev. A 54, 3948 (1996). 

V. A. Dzuba, M. G. Kozlov, S. G. Porsev, and V. V. 

Flambaum, Sov. Phys.-JETP 87, 885 (1998). 

V. N. Novikov and I. B. Khriplovich, JETP Lett. 22, 74 

(1975). 

N. H. Edwards, S. J. Phipp, P. E. G. Baird, and 
S. Nakayama, Phys. Rev. Lett. 74, 2654 (1995). 
V. V. Flambaum and I. B. Khriplovich, Zh. Exper. Teor. 
Fiz. 89, 1505 (1985). 



[13 

[i4; 
[is; 

[16 

[17; 
[is; 

[19 

[20 

[21 

[22 

[23 
[24 



C. Bouchiat and C. A. Piketty, Phys. Lett. B 269, 195 
(1991), Erratum:274,526 (1992). 

V. V. Flambaum and D. W. Murray, Phys. Rev. C 56, 
1641 (1997). 

W. C. Haxton, C.-P. Liu, and M. J. Ramsay- Musolf, 
Phys. Rev. Lett. 86, 5247 (2001). 

W. C. Haxton and C. E. Wie man, Ann. Rev. N ucl. Part. 

Sci. 51, 261 (2001), E-print: |nucl-th/0104026 . 

V. A. Dzuba, V. V. Flambaum, P. G. Silvestrov, and 

0. P. Sushkov, J. Phys. B 20, 3297 (1987). 

M. G. Kozlov, S. G. Porsev, and W. R. Johnson, Phys. 
Rev. A 64, 052107 (2001). 

1. B. Khriplovich, Parity non-conservation in atomic phe- 
nomena (Gordon and Breach, New York, 1991). 

S. G. Porsev, Y. G. Rakhlina, and M. G. Kozlov, Phys. 
Rev. A 60, 2781 (1999). 

S. G. Porsev and M. G. Kozlov, Phys. Rev. A 64, 064101 
(2001). 

V. G. Gorshkov, V. F. Ezhov, M. G. Kozlov, and A. I. 
Mikhailov, Sov. J. Nuc. Phys. 48, 867 (1988). 

D. E. Groom et al, Eur. Phys. J. C 15, 1 (2000). 
Note, that radiative corrections can change K2 quite no- 
ticeably. 



